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The order parameter of a quantum-coherent many-body system can include a phase degree of 
freedom, which, in the presence of an electromagnetic held, depends on the choice of gauge. Because 
of the relationship between the phase gradient and the velocity, time-of-Hight measurements reveal 
this gradient. Here, we make such measurements using initially trapped Bose-Einstein condensates 
(BECs) subject to an artificial magnetic held. Vortices are nucleated in the BEC for artihcial held 
strengths above a critical value, which represents a structural phase transition. By comparing to 
superhuid-hydrodynamic and Gross-Pitaevskii calculations, we conhrmed that the transition from 
the vortex-free state gives rise to a shear in the released BEC’s spatial distribution, representing a 
macroscopic method to measure this transition, distinct from direct measurements of vortex entry. 

Shear is also affected by an artihcial electric held accompanying the artihcial magnetic held turn-off, 
which depends on the details of the physical mechanism creating the artihcial helds, and implies 
a natural choice of gauge. Measurements of this kind offer opportunities for studying phase in 
less-well-understood quantum gas systems. 


While gauge invariance is central to our description of 
nature, specific physical situations lend themselves to a 
natural choice of gauge that matters. A Bose-Einstein 
condensate’s (BEC’s) order parameter includes a gauge- 
dependent phase. When a BEC is subjected to suffi¬ 
ciently rapid rotation (or a sufficiently strong artificial 
magnetic field [T] B), the BEC exhibits vortices: points 
at which the density vanishes and the phase is singu¬ 
lar. For a finite system, the structural phase transition 
from a state in which the phase varies smoothly to one 
with a single phase singularity occurs at a critical mag¬ 
netic field B CI , dependent on both particle-particle inter¬ 
actions and geometry [5|. In quantum gas experiments, 
images of time-of-flight (TOF) expanded clouds can re¬ 
veal the vortex-nucleation transition Em via the appear¬ 
ance of local minima in the imaged atomic density, each 
associated with a vortex core. Here, we detected this 
structural phase transition via an abrupt shape-change 
in TOF-expanded BECs accompanying the appearance 
of vortices. As we see below, this change depends upon 
the “natural” gauge choice for our experiment. 

In the description of many-body systems subject to 
a (real or synthetic) uniform magnetic field, the vector 
potential is generally written in the gauge most con¬ 
venient for the problem at hand. Although the sym¬ 
metric and Landau gauges are common choices, there 
is no a priori “natural” gauge. However, looking be¬ 
yond the simple specification that the magnetic field is 
uniform, the physical mechanism that creates the mag¬ 
netic field often suggests a specific gauge choice: any 
magnetic field created by a single electric current has a 
natural gauge determined by the geometry in which the 
vector potential is proportional to the electrical current. 
In this sense, the symmetric gauge is a natural gauge 
for an infinite cylindrical solenoid [Fig. 1(a)], and the 


Landau gauge is a natural choice for two parallel coun¬ 
terflowing sheets of current [Fig 1(b)]. Analogous to an 
infinite solenoid, the synthetic magnetic field for rotat¬ 
ing systems is most naturally expressed in the symmet¬ 
ric gauge A = (—Bye x + Bxe y )/2 because changes to 
the rotation rate lead to changes of A having the form of 
the symmetric gauge. As in the case for two sheet cur¬ 
rents, our experiment’s 0 synthetic magnetic field B is 
derived from an engineered vector potential [7] that most 
naturally takes the Landau form A = —Bye x , because 
the coupling between laser fields along e x only affects 
the component of A along e^, [Fig. 1(c)]. Two exper¬ 
iments starting with the same B and ending with zero 
magnetic field can manifest different physical outcomes 
without violating gauge invariance, since the electric field 
£ = — dA/dt associated with the change [2i is determined 
by the physical mechanism that creates B. In any gauge, 
the vector potential change A A resulting from modifica¬ 
tions to experimental parameters is, in this interpretation 
of the natural gauge, directly proportional to the natural 
gauge. 

In this work, we studied 87 Rb BECs subject to uniform 
laser-induced synthetic magnetic fields. The constituent 
bosons experienced a Lorentz force, just as would charged 
particles in a magnetic field. This synthetic magnetic 
field was continuously tunable from 0 to above B cx . BECs 
are characterized by the complex-valued order parame¬ 
ter ip( r) = y/ p( r) exp [i(p{ r)] with superfluid density p(r) 
and phase </>(r). The Gross-Pitaevskii equation (GPE) 
gives the time dependence of ip. The single-valuedness of 
the phase lends the superfluid system one of its defining 
properties: irrotationality. The current density J(r) = 
Me[^*(r)v^(r)], where v = [— ihV — A(v)\/m , describes 
the flow of particles and links the local velocity to its 
phase cp( r) via v(r) = [IiV(p( r) — .A(r)] /m for particles 
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(a) Solenoid: symmetric gauge 
(Rotating geometry) 



(b) Current sheets: Landau gauge 
(Raman geometry) 



(c) Experimental 
Raman geometry 



(d) 87 Rb level diagram 
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FIG. 1. (a) Solenoid and (b) counterflowing current-sheet geometries, shown with electrical current flow (red), vector potential 
A (blue), and magnetic field B (green). For each case, the natural gauge is illustrated: symmetric (suitable for rotating 
systems) A = (—Bye x + Bxe y )/2 in (a) and Landau (suitable for our Raman system) A = —Bye x in (b). (c) Experimental 
geometry. The elongated 8 ' Rb BEC (black) experienced an artificial field B oc e, created using counterpropagating Raman 
lasers (red and blue arrows), (d) Level diagram of the / = 1 ground state and excited 5 P states, showing two-photon Raman 
transitions. 


of mass m. Notice that irrotationality applies to the lo¬ 
cal per-particle canonical momentum p(r) = h'V(f>( r), 
i.e., V x p = 0, but neither to velocity nor mechanical 
momentum p m (r) = mv(r) = p(r) — .4(r). 

We prepared trapped, equilibrated [H] BECs in one of 
three configurations: (i) B = 0, a reference case described 
by standard superfluid hydrodynamics [Fig. 2(a,d)]; (ii) 
0 < B < B CI , described by modified superfluid hydro¬ 
dynamics [Fig. 2(b,e)]; and (Hi) B > B„, with vortices 
[Fig. 2(c,f)]. For all cases, we initiated TOF by abruptly 
removing the confining potential V (r) and rapidly mak¬ 
ing the vector potential ,4 = 0. Effectively [10], this 
rapid turn-off left cf>( r) unaltered and mapped the gauge- 
dependent canonical momentum p just before TOF (at 
t 0 - ) onto the gauge-independent mechanical momentum 
just after TOF began (at t 0 +): p(f 0 -) = Pm(^o+)- 

Measuring the shearing motion in TOF, which has 
contributions from the initial velocity field v(r) and the 
position-dependent electric force, allows us to distinguish 
between cases (i), (ii), and (in). 

Case (i): The B = 0 expansion of a repulsively interact¬ 
ing BEC released from a harmonic trap is well-studied: 
zero-point energy is typically negligible and interactions 
between the atoms dominate. In the ground state, the 
constant (in situ) phase gives p m (io+) = 0- The familiar 
inverted-parabola Thomas-Fermi profile of harmonically 
trapped BECs m is preserved by the interaction-driven 
expansion during TOF i~H? . 

Case (ii): Modest magnetic fields (0 < B < B cr ) al¬ 
ter the interaction-dominated TOF expansion, as was 
demonstrated in vortex-free rotating BEC systems [13]. 
The in situ phase <j>( r) is everywhere well-defined and 
the canonical momentum p(r) differs from the mechan¬ 
ical momentum p m (r) due to the presence of the vector 
potential that gives B = Vx A(r). The t = 0 canonical 


momentum in the Landau gauge 

P” = (y e * + xe v) W 

also defines p m (io+)> where —1 < e < lisa trap- 
and cyclotron-frequency-dependent anisotropy parame¬ 
ter [see Eq. (|5j) and Ref. m. Here the anisotropy is large, 
so that e is nearly — 1; compared to previous measure¬ 
ments [Tj in a cylindrically symmeteric trap (e = 0), the 
canonical momentum components are small. (In the cor¬ 
responding symmetric gauge expression e + 1 is replaced 

by e.) 

Case (Hi): Vortices significantly affect the TOF expan¬ 
sion. At the location of each vortex, the phase is singular 
[Fig. 2(c)], with a 27r-winding around it. As an example, 
a vortex centered at the origin of a cylindrically symmet¬ 
ric system contributes a phase <fi(r) = ta,nT 1 (y/x), and 
the t = 0~ canonical momentum in the Landau gauge is 
approximately 

p ... = ^Jpt- B (i±l) {ye , + xe ,y (2) 

As before, this in-trap canonical momentum becomes the 
mechanical momentum p m (io+) as TOF begins. 

We estimate the relative shear of these configurations 
from their momenta at t 0 +. Case (i) is shear free. 
For case (ii), the momentum along e x has typical scale 
|p m • e^l ss (e + l)By/2. In the presence of many vor¬ 
tices [case (Hi), in the diffused vorticity limit m] this 
increases by HN v y/R 2 , where R is the characteristic sys¬ 
tem size and N v is the total number of vortices. The 
difference is proportional to hn v y , where n v is the areal 
vortex density, leading to an abrupt increase in the shear¬ 
ing momentum when vortices enter [Fig. H-i)[. Mi¬ 
croscopically, this increased velocity originates from the 
spatial variations in phase associated with the vortices 
[Fig. 2(c)], 
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FIG. 2. GPE calculations of (a-f) in situ and (g-i) time- 
of-flight distributions. (a,b,c) Landau-gauge-m-sitw phase, 
(d,e,f) in situ density, and (g,h,i) TOF density profiles in the 
e^-ej, plane for B < B CI (left, centre) and B > B CI (right). 
The dashed grey ellipse in (a) indicates the Thomas-Fermi 
edge of the cloud; the phase is everywhere zero, (j-1) Exper¬ 
imental TOF images for systems subject to magnetic fields 
with cyclotron frequencies, flc = B/m. A small degree of 
shear is evident in (k), but the largest shear is found in (1), 
where the cloud is also fragmented, due to the vortices’ phase 
singularities. The dashed red lines in the lower images indi¬ 
cate the primary axis, to emphasize the small shear in (h) and 
(k). 


The transition between configurations ( ii ) and (Hi) 
at B cr occurs when the system can lower its en¬ 
ergy by admitting a vortex; this critical held depends 
on the system’s trap and interaction parameters [6]. 
The critical cyclotron frequency = B CI /m , above 
which vortices are energetically stable [HI dE; may 
be estimated as Og' lt = (5h/mR 2 ^)\n(R±/((). R± = 

[4/i/m(u; 2 + LUy)\ 7 is the mean transverse Thomas- 
Fermi radius in a harmonic potential with frequencies 
w x,y,z- The healing length £ = (fi 2 /2ra/z) 1/72 sets the 
characteristic vortex core size, where y, = a s /m) 

is the central mean-held energy, a s is the s-wave scatter¬ 
ing length El, and no is the central density. For the 
parameters in this experiment, VI c q/2tt = 12.0(1.1) Hz. 

We prepared BECs with N ss 1-4(3) x 10 5 in the 
/ = 1 ground state hyperhne manifold at the inter¬ 
section of two A = 1.064 fim laser beams [Fig. 1(c)]. 


The resulting potential was approximately harmonic 
and had measured frequencies {uj x ,(jJ v ,uj z }/ 2n = 

{10.1(1), 47.3(3), 90(1)} Hz [18]. We implemented an ar¬ 
tificial magnetic held H US] using the combination of 
two counterpropagating Ar = 790.1 nm Raman lasers 
[Fig. 1(d)] traveling along in conjunction with a 
(real) magnetic held B = ( B 0 + B'y)e y , giving a gradi¬ 
ent in detuning from Raman resonance gFBBB'/h along 
e y ranging from 0 to 640 Hz//zm. This gave effective 
cyclotron frequencies ranging from Qq/ 2ir = 0 to 20 Hz. 

To study the evolution of the density distribution af¬ 
ter mean-held-driven expansion, we released the atoms 
from the trap, and adiabatically transformed the Raman- 
dressed superposition into a single Zeeman level for imag¬ 
ing [ZD] in the hrst 2 ms of TOF [IB]. The cloud expanded 
for a total of 36.2 ms TOF before being imaged along e z 
[Fig.|:j,k,i)]. 

Figure [2ja-f ) shows computed in situ phase and den¬ 
sity distributions for a range of cyclotron frequencies. 
Vortices nucleate only above H cr , here only in panels (c) 
and (f). For this geometry, the system’s ground state 
consists of a linear chain of vortices; larger cyclotron 
frequencies or less anisotropy would result in a regular 
vortex lattice. (In the experiment, the vortices have not 
equilibrated to their ground state configuration i) Fig¬ 
ure [2](g-l) depicts calculated and measured TOF densi¬ 
ties; the cloud’s shear increases monotonically with in¬ 
creasing cyclotron frequency He = B/m [e.g., there is a 
small shear in Fig. ,k) and a large shear for Fig. IM)]- 
For strong artiheial helds the observed density distribu¬ 
tions were irregularly fragmented [Fig. [2jl)] , while the 
computed clouds were ordered. These observations are 
both consistent with the presence of vortices - disor¬ 
dered in the case of experiment - whose characteristic 
phase gradients cause density modulations after TOF. 
The generic ej,-aligned stripes present in TOF [Fig. §i>l)] 
result from the predominantly e y -expansion from the 
anisotropic trap. 

We obtained TOF images at various artificial field 
strengths, and fit the resulting 2D column densities to 
a sheared Thomas-Fermi profile 
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where A is the nJ OF density; R x , y are Thomas-Fermi 
radii; and a xy is a dimensionless shear parameter. We ob¬ 
tain a measure of the fragmentation from the fit residuals 
? = T,x,yl n 2 D as ( x ’y) - n 2 r,( x ,y)} 2 . Figure |Fa) shows 
that T sharply increases m around Hc/27t = 15 Hz, an 
indication of the transition to the BEC’s vortex phase. 

Additional evidence for the entrance of vortices into 
this anisotropically trapped BEC is the behavior of the 
shear parameter a xy : Fig. |3j[b) shows a xy sharply in¬ 
creasing above Hc/27t = 15 Hz, in reasonable agreement 
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0 5 10 15 20 

cyclotron frequency, Q. c /2n (Hz) 

FIG. 3. Signatures of vortex formation, (a) Fit residuals, T. 
indicating fragmentation due to vortices, (b) Shear parameter 
o xy . The experimental data points (grey circles) each reflect 
an average over tens of measurements, and the uncertainty 
bars denote the standard deviation of the mean [in (b) these 
are smaller than the symbol size]. Pink marks the region, in¬ 
cluding uncertainty, where the vortex nucleation transition is 
expected: Q.c/2tt = 12.0(1.1) Hz. The shear parameter cal¬ 
culated for irrotational BECs is denoted by the red dashed 
curve and the results of the Raman GPE calculations (dia¬ 
monds) connected by lines; neither has any free parameters. 
All uncertainties are statistical. 


with the predicted critical cyclotron frequency flg' lt /27r = 
12.0(1.1) Hz. We compare this result to two calculations. 

The first calculation, the usual hydrodynamic descrip¬ 
tion [153 , valid only in the absence of vortices, is 
modified to include the artificial magnetic field from a 
Landau-gauge vector potential [T5] , predicting the in situ 
density 



and the canonical momentum p“ [Eq. 01- Here, Ri = 
[2/z/md; 2 ] 1 ' 2 are modified Thomas-Fermi radii with effec¬ 
tive trapping frequencies O,, where <D 2 = w 2 + f2^[( 1 + 
e)/2] 2 , <D 2 = ujy + f2^[(l - e)/2] 2 , and ui z = uj z . The 
anisotropy parameter 


e 


~ 9 ~ 9 

U] x - UJ y 


w; 


- uit. 


( 5 ) 


can be obtained self-consistently from these equa¬ 
tions mi nu. These distributions are propagated in TOF 
using the hydrodynamic equations for comparison to the 


(a) Vector potentials (b) Calculated TOF signature 



position (pm) cyclotron frequency, Q. c /2% (Hz) 


FIG. 4. (a) Vector potentials used to calculate TOF response 
for cyclotron frequency Sic = 27 r x 16.7 Hz: the upper pan- 
nel shows the the vector potential corresponding to the exact 
Raman dispersion relationship (black) and the central-region 
fit to the Landau gauge (dashed blue) A x = —By (A y = 0 in 
both cases), while the lower panel shows the symmetric gauge 
components A x = —By/2 (blue dashed) and A y = Bx/2 (red 
dashed), (b) Shear parameter a xy computed via GPE for 
equivalent-held systems after TOF and H-turn-off for fields 
created by three physical mechanisms: Raman (blue); Lan¬ 
dau gauge (red); and symmetric gauge (green). 


measured distributions. From these calculated TOF dis¬ 
tributions, we extracted the shear parameter and com¬ 
pared it with the measurement, as shown in Fig. |3j[b) . 

The second calculation uses the GPE, which is valid 
for all our cases. We numerically found the GPE ground 
state at each B and then used the time-dependent GPE 
to calculate TOF evolution [18]. We find excellent agree¬ 
ment both with the amount of shear in the cloud for all 
configurations and in the location of the critical artificial 
field strength for vortex nucleation. 

To understand the impact of the experimentally dic¬ 
tated natural gauge, we investigated the effects of differ¬ 
ent gauge choices. We performed GPE calculations as 
described above to obtain a xy for two alternate gauge 
choices with the same B: the symmetric and Landau 
gauges. We compared these to the Raman system de¬ 
scribed above (which, as shown in Fig.[4](a), deviates only 
slightly from the ideal Landau gauge choice). Though the 
Landau and symmetric gauge atomic systems are identi¬ 
cal before release, their TOF responses differ because of 
their different electric field impulses at turn-off, which re¬ 
flect each system’s natural gauge [Fig.[4](b)]. For B < B CI , 
this electric field impulse produces a significantly larger 
shearing when the symmetric gauge is natural. This be¬ 
havior reflects the fact that the gauge-dependent in situ 
canonical momentum, which becomes p m (to+), is differ¬ 
ent in the two cases. Each system reacted to the entrance 
of vortices with a marked change in a xy , here positive, as 
set by the vortices’ direction of circulation, but with a 
significantly different character. 

In conclusion, we showed that the shape of the cloud 
observed after TOF is strongly affected by the presence of 
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vortices and signals the vortex nucleation transition. By 
exploiting the connection between phase and velocity, we 
were able to extract topological features of our system’s 
order parameter from the TOF density distribution. Our 
experiment’s anisotropic geometry and Landau-like nat¬ 
ural gauge led to a signal that made this transition par¬ 
ticularly clear. By comparing to what would have hap¬ 
pened in the symmetric gauge, we showed that the gauge 
choice can make a significant difference in the experi¬ 
mental outcome. In the future, techniques that match a 
system’s geometry to the artificial field’s natural gauge 
might be used to selectively excite edge modes in cold- 
atom quantum-Hall-like systems [2TH23] . 
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Experimental conditions 

Effects of modulation: This manuscript and Ref.ISlI 
share the same underlying data set, put to very different 
purpose. In this data set, the harmonic potential along 
e x was periodically modulated, but sufficiently weakly 
that the system remained in the linear response regime. 
In the data presented here, release always occurred at 
the same phase in the modulation cycle. (Modulation 
was used to determine transport coefficients, which evi¬ 
denced the superfluid Hall effect reported in Ref.ED) For 
the present analysis, we chose those modulation times for 
which the shear is entirely due to the induced circulation 
due to the applied TS and the vector-potential turn-off. 
We confirmed with our GPE simulations that the mod¬ 
ulation had no effect on the shapes of the clouds in situ 
or in TOF. 

Effects of Raman dressing: In addition to produc¬ 
ing an artificial magnetic field, the combination of Raman 
dressing and real magnetic field gradient used in these 
experiments contributed a scalar potential that locally 
depended on both Raman coupling strength and detun¬ 
ing from Raman resonance. (The spatially-dependent, 
lowest-energy eigenstates of the system with Raman cou¬ 
pling plus magnetic field and field gradient along e y ex¬ 
perienced an effective harmonic antitrapping potential.) 
This variation was measured experimentally by record¬ 
ing the frequency of dipole oscillation in the harmonic 
trap, confirming the weakened trapping potential along 
e y . We fit the resulting trap frequency to a phenomeno¬ 
logical second-order polynomial as a function of cyclotron 
frequency and included this modification of the trapping 
potential in the superfluid hydrodynamic calculations: 
ujy/2n = [c — a(fic/27r + &) 2 ], where a = 0.038(3) (Hz) -1 
and b = 2.1(9) Hz and c = 47.5(4) Hz . We also fully ac¬ 
counted for the Raman-induced anistropic effective mass 
in this system El E2] in the superfluid hydrodynam¬ 
ics calculations. In contrast, both the modified trap¬ 
ping potential and the effective mass are automatically 
accounted for in the Raman-system GPE calculations, 
where the full dispersion relationship was used, and no 
corrections are applied (see Notes on GPE calculations). 

Atom number: Our measurements of shear are rel¬ 
atively insensitive to atom-number differences between 
experimental realizations. For B < B CI , the shear a xy 
is depends only weakly on N. The atom number, N = 
1.4(3) x 10 5 was determined from the TOF Thomas-Fermi 
radii R y . The uncertainty is the standard deviation of the 
measurements, and represents shot-to-shot fluctuations. 


Trap and B turn-off 

The optical trapping potential was turned off at t = 0 
in less than 1 fj, s. Concurrently, the artificial mag¬ 
netic field was removed by adiabatically transforming 
the Raman-dressed superposition into a single Zeeman 
state, |/ = 1, mp = 1), by simultaneously ramping the 
Raman intensity to zero and sweeping the bias mag¬ 
netic field away from Raman resonance. This process 
was complete over the first 2 ms of TOF. Under this pro¬ 
cess, the 1/e time of turn-off of the vector potential was 
130 /is and during this time, the cloud expanded only 
slightly while artificial field effects were significant. Our 
GPE numerical calculations ignore the effects of this ar¬ 
tificial field during the initial mean-field expansion and 
assume that all expansion occurred without the Raman 
beams present. The agreement between our experimen¬ 
tal and numerical results validate this assumption, as ex¬ 
pected for an impulse that is short compared to typical 
timescales for expansion. In the main text, we simplify 
the turn-off and consider “t = 0 + ” to be the time when 
all Raman coupling effects are removed. 

Due to the large final Raman detuning, the final vec¬ 
tor potential Af = Af& x is non-zero, but uniform. This 
turn-off can be thought of as a two-step process. Af¬ 
ter the first step, A —> 0 and the mechanical momen¬ 
tum p^fo-i-) = p m (to_) — Bye x [equal to the initial 
canonical momentum p(fo)] was changed by an amount 
proportional to A in the Landau gauge. Physically, 
this change in mechanical momentum resulted from the 
electric field induced as A —> 0. In the second step 
A = 0 -+ A = A f , giving the final mechanical mo¬ 
mentum p m (to+) = Pm(^o+) — «4.f, shifted by a constant 
Af from the intermediate momentum p[ n (to+), because 
neither A = 0 nor A = Af have any spatial variation. 
As it is the spatially dependence that is of interest here, 
we consider in the main text the quantity p[ n (fo+)- Ef¬ 
fectively, the spatial variations of p^fo-e) and p m (f 0 +) 
are identical and independent of Af , allowing us to con¬ 
sider only the Af = 0 case in the main text, where we 
use p m (t 0 +). 

Notes on GPE calculations 

We modelled our system using a 2+1 dimensional sim¬ 
ulation, wherein we assumed a Thomas-Fermi profile 
along e 2 , and solved the resulting 2D GPE in the e x - 
e y plane, appropriately modified to account for the real 
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FIG. SI. GPE calculation of the time-of-flight distribution for 
a non-equilibrated vortex distribution with B > B CT ; besides 
equilibration time, all parameters are identical to those in 
Fig. 2(i). 
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3D profile by including a position-dependent interaction 
factor [S3]. We used imaginary-time propagation [54] 
to determine the initial wavefunctions, and evolved the 
system using a split-time spectral method jS5j- For the 
Raman-coupled GPE calculations, we assumed that the 
atoms remained in the lowest Raman dressed beam and 
used the energy versus momentum dispersion relation 
from the exact 3-level Raman-coupling Hamiltonian [S3] 
to describe the GPE’s kinetic energy (which in this case 
is only modified along e^,). In the presence of a detuning 
gradient along e y , this dispersion depended upon y. This 
simulation accounts for both the non-uniformity of 3 and 
all contributions to the potential energy. The calculation 
is valid for all three cases described in the main text, and 
for example correctly predicts the low-field shear, VLq for 
the onset of vortex nucleation, and the shear with vor¬ 
tices. 

The regular structure in the simulated TOF distribu¬ 
tion we see in case (Hi) [Fig. 2(i)] results from the reg¬ 
ularity of vortex positions in situ. If we start the GPE 
calculation with vortices seeded at random positions and 
allow the system to equilibrate, the vortices will even¬ 
tually form a single row along the long axis, as seen in 
Fig. 2(f). Figure SI shows a simulation with the same 
parameters as Fig. 2(i), but before equilibration. This 
strongly resembles the experimental result with the same 
parameters, Fig. 2(1). 

In our comparisons to the idealized gauge cases, we 
fit the vector potential computed for the Raman-coupled 
Hamiltonian near the trap centre to that of the Landau 
gauge for a uniform 3. We then used the same field value 


3 for the symmetric gauge calculation. 

Note that the cyclotron frequency differs by a factor of 
two from the rotation frequency used in a similar expres¬ 
sions for rotating systems, as in Ref. lS6l and elsewhere. 


Vortex nucleation 


In the main text, we approximate the canonical mo¬ 
mentum in case (Hi) under the assumption of a single 
vortex at the centre of the cloud, and make the assump¬ 
tion that the healing length, i.e., the vortex core size, is 
much smaller than the Thomas-Fermi radius. 

Fragmentation parameter: The measured fragmen¬ 
tation parameter plotted in Fig. 3(a) in the regime of case 
(Hi) is non-monotonic as a function of flc due to the de¬ 
tails of the imaging process: though we see an initial 
increase as B increases beyond B CI , it does not continue 
to increase. As more vortices enter the system, the TOF 
density variations increase in their spatial frequency, and 
due to the limited resolution of our imaging, this measure 
of fragmentation decreases for higher spatial frequencies. 
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